Atomic Fock states by gradual trap reduction: from sudden to adiabatic limits 
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We investigate the possibility to form high fidelity atomic Fock states by gradual reduction of 
a quasi one dimensional trap containing spin polarized fermions or strongly interacting bosons in 
the Tonk-Girardeau regime. Making the trap shallower and simultaneously squeezing it can lead to 
the preparation of an ideal atomic Fock state as one approaches either the sudden or the adiabatic 
limits. Nonetheless, the fidelity of the resulting state is shown to exhibit a non-monotonic behaviour 
with the time scale in which the trapping potential is changed. 
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I. INTRODUCTION 

Preparation of atomic states containing exactly a fixed 
number M of atoms (Fock states) is of importance for a 
wide range of applications, from studying ultracold chem- 
istry and few-body physics to precision measurements 
and quantum information processing. The aim of such a 
preparation is to create a quantum state with the mean 
number of atoms {n) equal to the chosen M and its mean 
variance as small as possible. Different proposal based on 
a time-dependent modulation of the confining potential 
have been put forward both for optical traps [IHS] and 
optical lattices [3 [7] . 

Recent experiments have shown the possibility of 
achieving atom-number sub-Poissonian statistics in a 
quantum degenerate gas with repulsive interactions [S]. 
The key idea is that the mean number of trapped atoms 
can be controlled by adiabatically reducing the depth of 
the trap while expelling the excess of atoms. The pre- 
cision of the control improves by maximising the energy 
splitting between states with different number of particles 
which ultimately allows one to discriminate such states 
by modulating the depth of the trapping potential. This 
technique is referred as to atom culling [1] . For ultracold 
gases confined in tight-waveguides, it works optimally in 
the strongly interacting regime for Bosonic samples 11- 
H], i. e., in the Tonks- Girardeau gas, where the repulsive 
interactions lead to an effective Pauli exclusion principle 
[S]. Alternatively, Fock state preparation is optimized 
as well with a spin-polarized non-interacting Fermi gas 
[S]. Bearing in mind that both systems (which are 
dual under Bose- Fermi mapping [9J) share all local cor- 
relation functions, and in particular the atom-number 
distribution [12] . we shall address the polarized fermions 
for brevity in the following. Nonetheless, all results in 
this paper apply to both systems. 

Non-interacting spinless fermions when placed in a 
trap occupy the lowest one-particle levels in the ground 
state. Changing the shape of the potential trap aims to 
expel redundant atoms into continuum levels leaving the 



bound states of the final trap filled to its maximum ca- 
pacity, (n) = M . This condition together with the Pauli 
exclusion principle ensures high fidelity of the prepara- 
tion: since no more than M atoms can be distributed over 
M levels of the final well, the variance = {'h?) — (n)^ 
must vanish, and precisely M atoms will be found in each 
individual case. A slow change of the trapping potential 
would guarantee, by virtue of the adiabatic theorem, full 
occupation of the states in the final well. The adiabatic- 
ity may, however, require times large compared to the 
time scales typically involved in ultracold atom experi- 
ments and one might wish for a faster way to achieve 
full final occupation. One of the counter intuitive re- 
sults obtained in Refs. |2l [3] is that an infinitely fast 
"sudden" change provides an alternative to the adiabatic 
route which can lead to the preparation of ideal Fock 
states, when making the well shallower (weakening of the 
trap) is accompanied by also making it narrower (squeez- 
ing). Further, at zero temperature, the best results are 
obtained 2j in case the projector on the subspace of all 
filled initial states Aq contains the projector A on the 
subspace spanned by all bound states of the final well 
(A C Aq), i.e., when on that subspace Ag can be ap- 
proximated by unity. In addition, non-zero temperature 
effects can be overcome by starting with a larger initial 
sample [3]. 

However, the sudden limit may not be accessible in a 
given experimental set up, and the main purpose of the 
present paper is to explore the behaviour of the trapped 
atoms when a change of the potential is neither sudden, 
nor adiabatically slow. In particular, we will show that 
the efficiency of the Fock state preparation exhibits a 
non-monotonic behaviour as a function of the duration 
of the quench of the trapping potential. We will also 
address other questions such as the dependence of final 
occupation on the depths of the initial and the final wells, 
and the escape time required for the expelled atoms to 
leave the well region. The rest of the paper is organised 
as follows. In Sect. II we will study the dependence of 
the final occupation on the time in which the trapping 
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potential changes its shape. In Sect. Ill we briefly dis- 
cuss the behaviour of the variance when the final well is 
filled close to its maximal capacity. In Sections IV we 
review two different mechanisms responsible for the for- 
mation of nearly ideal Fock states in the adiabatic and 
sudden limits. In Sections V we discuss the escape of 
expelled particles from the final well region, and Section 
VI contains our conclusions. 



II. ATOM COUNTING STATISTICS 
FOLLOWING A CHANGE OF THE TRAPPING 
POTENTIAL 

Consider a number of non-interacting spinless fermions 
(fermionised bosons) initially trapped in a one-dimesional 
rectangular square well of depth Vi and width Li. The 
potential V{x, t) undergoes, over a time T, a transforma- 
tion such that its final shape is also a rectangular well, 
but with new depth and width, V/ and Lf, respectively. 
Although other types of evolution are possible, we will 
consider linear change of both the depth and width of 
the trap, 

V{t) ^V, + [Vf - V,)t/T, (1) 
Lit)=L, + {Lf-L,)t/T, 

where the ramping time T determines whether evolution 
of the potential is rapid or slow. Our aim is to control 
the mean number {n{T)) of the fermions trapped in the 
final trap, while minimizing the variance 

<(T) ^ {n'iT)) (r^(^))^ (2) 

for the case of zero temperature. As already mentioned in 
the introduction, non-zero temperature effects turn out 
to be non-critical and can be conveniently overcome by 
starting with a larger initial sample as described in [3]. 
We will further assume that the initial and final traps 
support K and M < K bound states, respectively, and 
that there are N < K fermions occupying the first N 
levels of the initial well. A well known technique ^TQ\ - 
[T^ allows to express {n{T)) and (j]^{T) for a system of 
non-interacting fermions in terms of the solutions of the 
corresponding one-particle Schrodinger equation. Let us 
denote by (j>'^{x) {(p-^ (x)) the states of the initial (final) 
trap. Bound (scattering) states will be labeled by a Latin 
(Greek) letter as a subindex. In the dimensionless vari- 
ables x/Li and t/to, where 

to = fiL^/h, (3) 

and ^ is the atomic mass, the time-dependent single- 
particle eigenstates obey 

iMl, = -^dlc^:, + W{x,t)ct>l,, n=l,...,N, (4) 

with W{x, t) = V{x, t)tolh. For a trap size of 80 /zm the 
time to for Rb and Cs atoms takes values of 8.8 s and 



13.4 s, respectively. Following [51 [31 [TOl E] we obtain 

M 

(n(T))=^(0^^|AT|0,^) (5) 

and 

M 

a%{T) = (n(T)) - ^ (0^^|AtAAt)|</.J). (6) 

Here, A is the projector on the subspace spanned by the 
one particle bound states j — l,2,...,Af of the 

final well, 

M 

Similarly, A^ is the projector onto the subspace spanned 
by the orthogonal states obtained by the time evolution 
of the one particle bound states \4>n), n = 1,2, . . . , N, in 
the initial well, 

N 

AT^Y.\^n){€\, \cl^l)^U{TMl), (8) 
Ji=l 

with U (T) denoting the evolution operator corresponding 
to Eq. Q. Note that Eqs. ^ and ^ are generalisations 
of Eqs. (9) and (10) obtained in Ref. [2^ for the sud- 
den limit, with the initial bound states replaced by time 
evolved states We further notice that knowledge of 
time evolved states ([8| allows one to obtain the full atom- 
number distribution p{n) from the characteristic func- 
tion F{9) = Tr[/5e'^^"'^], as a Fourier transform, p{n) = 
^ C^e-'"'''FiS)de, with n = 1,...,M ^MM- Using 
the projector for the bound subspace in the final config- 
uration, F{e) = detAwith A = [l-|-(e*^-l)AAT]. In the 
basis of single-particle eigenstates \4>m), the elements of 
the matrix A read Anm = <5nm + [exp(i0) — l]((/)^|AT|(?!>m). 

With the problem reduced to numerical evaluation of 
the corresponding one particle states, we employ the 
Crank-Nicolson method to solve Eq. Q with the initial 
conditions (t)\^{x,t — 0) n = 1,2. ..N and zero boundary 
conditions at the edges of the numerical box, x = ±Lf,. 
To avoid reflections from the boundaries we introduce an 
absorbing potential proposed by Manolopoulos [TH [T3] 
(see appendix A for details). Following Refs. [21 [3], we 
will refer as "weakening" to the case where the trap is 
made shallower while keeping its width constant, Vi > 
Vf, Li = Lf , while reducing the trap's width will be 
called "squeezing". 

From Eqs. ([5| and ^ it is easy to deduce [51|3] that an 
ideal Fock state with (n(T)) = M and aj^{T) = would 
be prepared provided that A C Ay, this is, the space 
spanned by the state before the end of the quench should 
enclose that spanned by the Fock state to be prepared 
[IB] . In the following we shall describe different physical 
implementations to fulfill this requirement. 
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III. TWO WAYS TO ARRIVE AT A GOOD 
FOCK STATE 




FIG. 1. Dependence of the final occupation on the ramping 
time (pure weakening): {n(T)) (solid) and {n{T)} — aN{T) 
(dashed) vs r = T/U. The potentials are Wi = (19.5)^7r^/2, 
Wf = (4.95)^7r'^/2, corresponding to maximum capacities 
K = 2Q and M = 5, respectively, and JV = if = 20, it is, all 
initial levels occupied. 



In the adiabatic limit, where the change of the potential 
shape is slow, we can expect the first M bound states of 
the initial well to follow into the M bounds states of the 
final well 



\4>J,) = eMi<^nMi), n = 1,2,...,M, 



(9) 



where <&„ is a real phase. Thus, on the subspace of final 
bound states we have 



M 
n=l 



M 



(10) 



Figure [T] shows the final occupation {n{T)) and its vari- 
ance as functions of the ramping time for the case of pure 
weakening of the trap. As expected, the final occupation 
increases and the fidelity improves as one approaches the 
adiabatic limit, where the variance of the resulting state 
vanishes. This result (a consequence of the adiabatic the- 
orem) does not depend on whether weakening, squeezing 
or a combination of both techniques is applied. 

It is worth noting that, to our knowledge, the question 
of exactly how slowly the potential must evolve to ensure 
adiabaticity in the case the M-th level in the final well lies 
close to the edge of continuum remains open and requires 
further investigation. As mentioned in the introduction, 
the applicability of the adiabatic method may be limited 
by a finite lifetime of the trapped condensate. 

In the sudden limit where the shape of the well changes 
almost instantaneously we have U{T) « 1, = {(pl) 



and = i. e. At = Aq. Let us recall that 

if A C At, i.e, the final well would contain a Fock state 
with precisely M fermions. Remarkably, the use of either 
either weakening and squeezing in the sudden limit lead 
to quasi-Fock states of limited fidelity [2] , while the com- 
bination of both techniques can lead to ideal Fock states 
[21 E] . Generally, sudden quenches of the trap might lead 
to undesirable excitations of the transverse modes break- 
down the effective one-dimensional character of the sys- 
tem. 



IV. GRADUAL MODULATION OF THE TRAP 
POTENTIAL 

We have seen that both adiabatic and sudden changes 
of the trapping potential can in principle lead to the 
preparation of ideal Fock states, but as idealized limits 
might be of limited relevance to experimental implemen- 
tation of atom culling techniques. 

Motivated by this observation we next look at the mean 
particle number and variance of a state resulting from a 
combined process of squeezing and weakening of the trap 
in finite time, see Fig. [2] Here, (n(T)) stays close to 
the maximum capacity of the well for all ramping times, 
but exhibits a dip at T/to « 0.03 — 0.1. Notice that the 
adiabatic preparation relies on the adiabatic following of 
the first M states of the initial trap onto the M states 
of the bound subspace of the final trap T]. By contrast, 
the sudden preparation relies on the instantaneous reso- 
lution of the target state \M) within the space spanned 
by the initial state [21 [3] . The non-monotonic behaviour 
of (?^(^)) reported in Fig. |2] results from a simultaneous 
failure of these two different mechanisms that allow the 
perfect resolution of the desired Fock state in either of 
the T — >■ 0, oo limits. 

We further note (see Fig. 2, dashed) that in the case 
of combined weakening and squeezing the mean number 
of atoms and its variance (both dimensionless, of course) 
add almost exactly to the maximum capacity of the final 
well, M, for all ramping times, and in the next section 
we discuss this approximate "sum rule" in some detail. 



V. SUB-POISSONIAN STATISTICS AT 
NEARLY FULL FINAL OCCUPATION 

Let us start by introducing the deficiency operator 

i = A-AT. (11) 
Inserting (111 into ^ and ([6| yields 



M 



.(r)) = Af-^(</)J|i|<^^^) 



(12) 



and 



M 



M 



aUT) = J2(^i\^\^j)- E \i<f'i\M)\'- (13) 

j = l j,k = l 



4 



4,96 



C 4,92 
V 



4,88 

4,84 
0,02 

0,01 



merits of the distribution p, 





0,02 




1 


j\x10^ 


0.015 

2^ 0,01 
Q. 

0.005 


I =0.005 


1=0.095 















0,05 


0,1 0,15 0,2 
T 



0,1 



0,2 0,3 
T 



0,4 



0,5 



FIG. 2. a). Same as Fig. [Tlbut for weakening combined with 
squeezing, Lf/Li — 0.6. ^Iso shown by a dot-dashed line is 
the l.h.s. of Eq. flS] ). 

b) Probabilities p(4) (solid) and p(3) x 10^ (dashed) in Eqs. 
(16 1 vs T = T/to- The inset shows the probability p{4) and 
the ramping time interval for which p(4) > 0.5%. 



From Eqs. (12) and (13) follows a "sum rule" 



M 

{n{T))+al{T)^M- ^ (14) 

Note that the last term is quadratic in the operator A — 
1/ — At so that if A is "small" [16] the variance can be 
expressed in terms of {n(T)) through the approximate 
relation 



{n{T))+aUT)^M, 



(15) 



as demonstrated in Fig. 2a. This is an example of non- 
Poissonian behaviour resulting from indistinguishability 
of the fermions involved. Indeed, for a Poissonian distri- 
buti on o ne would expect ct^(T) — {n{T)) whereas from 
Eq. ^ we have cr^(T) = M- (n(T)) < {n{T)). 

Note also that in the case of almost full final occu- 
pation, knowing {n(T)) and cr^(T) allows one to recon- 
struct the full counting statistics. Thus, neglecting the 
probabilities to trap less than M — 2 fermions, p{k) « 
for = 0, 1, . . . ,M-3, we can obtain p(Af- 2), p{M-l) 
and p{M) from the normalisation and the first two mo- 



M 



m = 0,1,2. 



(16) 



k=M~2 



Results of such a reconstruction are shown in Fig. 2b. 
For the set corresponding to Fig. [2}3 the probability to 
trap three fermions, p(3), is negligibly small, and one 
would obtain 4 fermions at most in 1.2% of all cases for 
T = T/io ~ 0.01. To reduce this fraction to 0.5% one has 
an option to choose the ramping time either less than 
Ti < 0.005 or longer than T2 > 0.095 shown on the inset 
in Fig. [2]d. 



VI. ESCAPE OF EXPELLED PARTICLES 
FROM THE FINAL TRAP REGION 



Although the final occupation of the bound states in 
the final well is determined by the time T after which 
the potential no longer changes, the obtained Fock state 
can only be used after the atoms expelled into the con- 
tinuum have left the well area. To obtain an estimate for 
how fast this would happen we have chosen an interval 
51: —l.^Lf < X < 1.5Lf, containing the final well region, 
and monitored the mean number of atoms in Q, (no(t)), 
during the ramping, t < T, as well as for t > T. Figure 
|3] shows the results for combined weakening and squeez- 
ing, Wf/W, = 0.18, Lf/Li = 0.6. for T/to = 0.05 and 
T/to = 1-0, close to the sudden and the adiabatic limits, 
respectively. With this choice of parameters, the initial 
well contains 20 atoms and the final well's maximum ca- 
pacity is 5. 

In the nearly sudden limit, T/to — 0.05, a considerable 
fraction of expelled atoms remain in the region by the 
time the well achieves its final configuration. One has 
then to wait a duration of the order of to for the mean 
number of atoms in to settle to its final value {n{T)) ss 
M = 5. A more detailed analysis shows that it is the 
M + 1 (in this case, the sixth) state of the initial well 
which is delayed most in leaving the area. We note also 
that this decay is not exponential and, therefore, cannot 
be attributed solely to trapping of an atom in one of the 
resonances supported between the edges of the final well. 

In the nearly adiabatic case, T/to = 1; initial one- 
particle levels are gradually pushed into the continuum, 
and most of the atoms have time to leave while the 
potential shape is still changing. Once the change has 
stopped one still has to wait approximately ig for the 
contribution from the M + 1-th state to clear the area. 
Although comparable in magnitude, the wait is some- 
what longer than that in the nearly sudden limit, mostly 
because the expelled atoms receive more energy and move 
faster if the change of the potential shape is sudden. 
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Appendix A: Complex absorbing potential 

The absorbing potential employed in this paper has 
the form suggested in [T31 [H]. With the wavefunction 
required to vanish at the edges of the computational box, 
X = ±Lft denoting the size, the absorbing potential is 
chosen to be zero for < x < Lf,/2. For Lf,/2 < x < Lf, 
it is given by 



Vabs (x) 



-iD 



Az - Bz-^ 



{C-zf {C + zf 



(Al) 



where z = C{2x/U - 1), C = 2.62206, A = (1 - IG/C^), 
B = (1 - 17/C3)/C2, and D = C'^Ll/1.28. Finally, 
Vabs{x) — Vabsi—x). The value L^/Li — 10 was used 
and the Schrodinger equation was solved numerically on 
a grid of 4 • 10^ points spanning the interval [-~Li,,Lb]. 



In summary, we have investigated formation of atomic 
Fock states by gradually changing the trapping potential 
and expelling the excess of atoms out of the trap. For 
the case in which either the depth (weakening) or the 
width (squeezing) of the potential trap are reduced, an 
increase of the ramping time leads to an improvement in 
the fidelity of the procedure until the adiabatic limit is 
reached. This monotonic behaviour in a relevant vari- 
able is generic in phenomena with sudden to adiabatic 
crossovers. For example, in a Landau-Zener transition in 
which the crossing time is gradually increased, the final 
state goes from one level to the other depending on that 
time and the corresponding degree of adiabaticity of the 
passage. However, in our case, a combination of both 
procedures (squeezing and weakening), which produces 
high fidelity Fock states in the sudden limit, remains su- 
perior to either of both operations for all ramping times. 
This opens a new route for feasible atomic Fock state cre- 
ation by trap reduction, approximating the sudden limit 
in the laboratory. The behaviour along the crossover is 
non monotonic and for finite quenching times one ob- 
serves a decrease in the mean number of trapped atoms 
which is accompanied by a corresponding increase in its 
variance. Therefore in this combined scheme we find the 
rather unusual property that both sudden and adiabatic 
limits lead to the same results, certainly through differ- 
ent mechanisms, with a non-montonic decrease of fidelity 
in between. 



Appendix B: Final occupation vs the depth of the 
initial trap 

The discussion in Section IV suggests that, close to 
the sudden limit and for weakening and squeezing ap- 
plied together, starting with a deeper well fully filled, 
N = K, should improve fidelity. In addition, starting 
with a deeper initial well with a fixed number of level 
filled, N, would not be as beneficial, since the frequency 
of oscillations of the first unfilled state, (x|05y+i): is al- 
most independent of Vi. This is illustrated in Fig. |4] 
where the solid line shows the final occupation of the 
well with the maximum capacity M — 5 with all the lev- 
els of the initial well filled, N — K. The graph show 
peaks which correlate with initial depths at which a new 
level appears in the initial well. In the same figure, the 
dashed line shows the dependence of {n{T)) on Vi when 
only first 25 levels of the initial well are occupied. In this 
case, (n(r)) is no longer sensitive to the appearance of 
new bound states, but yields a slightly lower fidelity of 
the preparation. 

The case of pure weakening shown in Fig. [5] is more 
complex. There increasing the depth of a fully filled well 
leads to a decrease in {n{T)) (solid), which is also highly 
sensitive to the increase in the number of initial bound 
states (see inset in Fig. [s]). As in Fig. |4j partial filling of 
the initial well removes the structure in the dependence 
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FIG. 4. Final occupation in the sudden limit, {n{T — 0)), 
for the case of weakening with squeezing vs. the dimen- 
sionless initial well depth for N = K (solid) and N — 25 
(dashed). Other parameters are Wf = (4.95)^7r^/2(L//Li)^ 
(M — 5) and Lf/Li = 0.6. Note that for a deep well 
(2Wy^^/n gives an estimate of the number of bound states, 
K « Int{{2Wy/^/Tv). 
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FIG. 5. Same as Fig. |4]but for the case of pure weakening, 
Li = Lf and ior N ^ K (solid) and = 25 (dashed). The 
insert shows the correlation between the structure in the A'' = 
K curve and the number of the bound states supported by 
the initial well. 



of {n{T)) on Vi but leads to lower final occupation in the 
deep well limit. 
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